## Loading required package: foreach
##
## Attaching package: 'foreach'
## The following objects are masked from 'package:purrr':
##
## accumulate, when
## Loading required package: iterators
##
## Attaching package: 'data.table'
## The following objects are masked from 'package:lubridate':
##
## hour, isoweek, mday, minute, month, quarter, second, wday, week,
## yday, year
## The following objects are masked from 'package:dplyr':
##
## between, first, last
## The following object is masked from 'package:purrr':
##
## transpose
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`
## [1] 1

## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.060185 -0.003642 0.000510 0.003365 0.041035
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -0.0022339 0.0001474 -15.16 <2e-16 ***
## test_set$FCH4 0.1299140 0.0073970 17.56 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.006975 on 2291 degrees of freedom
## Multiple R-squared: 0.1187, Adjusted R-squared: 0.1183
## F-statistic: 308.5 on 1 and 2291 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.170771 -0.002274 0.000127 0.001837 0.178992
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -1.229e-03 6.965e-05 -17.65 <2e-16 ***
## result$FCH4 5.134e-01 3.089e-03 166.22 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.006624 on 9169 degrees of freedom
## (18521 observations deleted due to missingness)
## Multiple R-squared: 0.7508, Adjusted R-squared: 0.7508
## F-statistic: 2.763e+04 on 1 and 9169 DF, p-value: < 2.2e-16
##
## [1] 1
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] 2

## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.071636 -0.003196 0.000791 0.003452 0.103163
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -0.0026030 0.0001583 -16.44 <2e-16 ***
## test_set$FCH4 0.1006344 0.0071447 14.09 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.00754 on 2291 degrees of freedom
## Multiple R-squared: 0.07969, Adjusted R-squared: 0.07929
## F-statistic: 198.4 on 1 and 2291 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.194724 -0.002160 0.000432 0.002121 0.085284
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -1.606e-03 7.047e-05 -22.79 <2e-16 ***
## result$FCH4 4.566e-01 3.125e-03 146.12 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.006702 on 9169 degrees of freedom
## (18521 observations deleted due to missingness)
## Multiple R-squared: 0.6996, Adjusted R-squared: 0.6995
## F-statistic: 2.135e+04 on 1 and 9169 DF, p-value: < 2.2e-16
##
## [1] 2
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] 3

## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.040536 -0.003163 0.000527 0.002967 0.045172
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -0.002112 0.000138 -15.31 <2e-16 ***
## test_set$FCH4 0.101434 0.005867 17.29 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.006558 on 2291 degrees of freedom
## Multiple R-squared: 0.1154, Adjusted R-squared: 0.115
## F-statistic: 298.9 on 1 and 2291 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.161699 -0.002215 0.000296 0.001926 0.090340
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -1.402e-03 6.822e-05 -20.56 <2e-16 ***
## result$FCH4 4.409e-01 3.026e-03 145.71 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.006488 on 9169 degrees of freedom
## (18521 observations deleted due to missingness)
## Multiple R-squared: 0.6984, Adjusted R-squared: 0.6984
## F-statistic: 2.123e+04 on 1 and 9169 DF, p-value: < 2.2e-16
##
## [1] 3
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] 4

## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.074293 -0.003247 0.000909 0.003260 0.078962
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -0.0028446 0.0001404 -20.27 <2e-16 ***
## test_set$FCH4 0.1094572 0.0053669 20.39 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.0067 on 2291 degrees of freedom
## Multiple R-squared: 0.1537, Adjusted R-squared: 0.1533
## F-statistic: 416 on 1 and 2291 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.168504 -0.002251 0.000502 0.002290 0.075533
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -1.861e-03 6.816e-05 -27.3 <2e-16 ***
## result$FCH4 4.018e-01 3.023e-03 132.9 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.006482 on 9169 degrees of freedom
## (18521 observations deleted due to missingness)
## Multiple R-squared: 0.6584, Adjusted R-squared: 0.6584
## F-statistic: 1.767e+04 on 1 and 9169 DF, p-value: < 2.2e-16
##
## [1] 4
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] 5

## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.055604 -0.003821 0.000448 0.003576 0.054082
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -0.0025380 0.0001681 -15.09 <2e-16 ***
## test_set$FCH4 0.1029926 0.0085694 12.02 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.007965 on 2291 degrees of freedom
## Multiple R-squared: 0.05931, Adjusted R-squared: 0.0589
## F-statistic: 144.4 on 1 and 2291 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.191550 -0.002307 0.000214 0.001947 0.194081
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -1.309e-03 7.525e-05 -17.39 <2e-16 ***
## result$FCH4 5.135e-01 3.337e-03 153.87 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.007157 on 9169 degrees of freedom
## (18521 observations deleted due to missingness)
## Multiple R-squared: 0.7208, Adjusted R-squared: 0.7208
## F-statistic: 2.368e+04 on 1 and 9169 DF, p-value: < 2.2e-16
##
## [1] 5
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] 6

## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.059442 -0.003249 0.000606 0.003189 0.045591
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -0.002504 0.000144 -17.38 <2e-16 ***
## test_set$FCH4 0.107270 0.006455 16.62 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.006867 on 2291 degrees of freedom
## Multiple R-squared: 0.1076, Adjusted R-squared: 0.1072
## F-statistic: 276.2 on 1 and 2291 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.111727 -0.002138 0.000367 0.002146 0.091988
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -1.635e-03 6.698e-05 -24.41 <2e-16 ***
## result$FCH4 4.512e-01 2.971e-03 151.88 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.00637 on 9169 degrees of freedom
## (18521 observations deleted due to missingness)
## Multiple R-squared: 0.7156, Adjusted R-squared: 0.7155
## F-statistic: 2.307e+04 on 1 and 9169 DF, p-value: < 2.2e-16
##
## [1] 6
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] 7

## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.066273 -0.003303 0.000721 0.003176 0.037723
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -0.0023071 0.0001374 -16.79 <2e-16 ***
## test_set$FCH4 0.1308690 0.0060818 21.52 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.006534 on 2291 degrees of freedom
## Multiple R-squared: 0.1681, Adjusted R-squared: 0.1678
## F-statistic: 463 on 1 and 2291 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.176707 -0.002282 0.000344 0.002076 0.095323
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -1.485e-03 6.467e-05 -22.97 <2e-16 ***
## result$FCH4 4.479e-01 2.868e-03 156.18 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.00615 on 9169 degrees of freedom
## (18521 observations deleted due to missingness)
## Multiple R-squared: 0.7268, Adjusted R-squared: 0.7268
## F-statistic: 2.439e+04 on 1 and 9169 DF, p-value: < 2.2e-16
##
## [1] 7
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] 8

## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.045952 -0.003452 0.000762 0.003332 0.061767
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -0.0024565 0.0001495 -16.43 <2e-16 ***
## test_set$FCH4 0.1202037 0.0060581 19.84 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.007116 on 2291 degrees of freedom
## Multiple R-squared: 0.1466, Adjusted R-squared: 0.1463
## F-statistic: 393.7 on 1 and 2291 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.147316 -0.002235 0.000376 0.002097 0.085438
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -1.544e-03 6.936e-05 -22.26 <2e-16 ***
## result$FCH4 4.350e-01 3.076e-03 141.40 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.006596 on 9169 degrees of freedom
## (18521 observations deleted due to missingness)
## Multiple R-squared: 0.6856, Adjusted R-squared: 0.6856
## F-statistic: 1.999e+04 on 1 and 9169 DF, p-value: < 2.2e-16
##
## [1] 8
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] 9

## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.043495 -0.003283 0.000498 0.003219 0.068228
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -0.0026574 0.0001369 -19.4 <2e-16 ***
## test_set$FCH4 0.0952281 0.0052045 18.3 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.006536 on 2291 degrees of freedom
## Multiple R-squared: 0.1275, Adjusted R-squared: 0.1271
## F-statistic: 334.8 on 1 and 2291 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.168388 -0.002337 0.000424 0.002329 0.080888
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -0.0017827 0.0000703 -25.36 <2e-16 ***
## result$FCH4 0.4039347 0.0031179 129.55 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.006686 on 9169 degrees of freedom
## (18521 observations deleted due to missingness)
## Multiple R-squared: 0.6467, Adjusted R-squared: 0.6467
## F-statistic: 1.678e+04 on 1 and 9169 DF, p-value: < 2.2e-16
##
## [1] 9
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] 10

## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.048676 -0.003339 0.000558 0.003213 0.108503
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -0.0023348 0.0001481 -15.77 <2e-16 ***
## test_set$FCH4 0.1153319 0.0070928 16.26 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.007024 on 2291 degrees of freedom
## Multiple R-squared: 0.1035, Adjusted R-squared: 0.1031
## F-statistic: 264.4 on 1 and 2291 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.154339 -0.002194 0.000272 0.001979 0.096263
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -1.392e-03 6.634e-05 -20.99 <2e-16 ***
## result$FCH4 4.698e-01 2.942e-03 159.67 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.006309 on 9169 degrees of freedom
## (18521 observations deleted due to missingness)
## Multiple R-squared: 0.7355, Adjusted R-squared: 0.7355
## F-statistic: 2.55e+04 on 1 and 9169 DF, p-value: < 2.2e-16
##
## [1] 10
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] 11

## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.058739 -0.003603 0.000674 0.003438 0.089687
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -0.0026711 0.0001555 -17.18 <2e-16 ***
## test_set$FCH4 0.1062575 0.0083367 12.75 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.007369 on 2291 degrees of freedom
## Multiple R-squared: 0.06621, Adjusted R-squared: 0.06581
## F-statistic: 162.5 on 1 and 2291 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.160146 -0.002130 0.000347 0.001977 0.129836
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -1.455e-03 6.684e-05 -21.77 <2e-16 ***
## result$FCH4 4.932e-01 2.964e-03 166.37 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.006357 on 9169 degrees of freedom
## (18521 observations deleted due to missingness)
## Multiple R-squared: 0.7512, Adjusted R-squared: 0.7511
## F-statistic: 2.768e+04 on 1 and 9169 DF, p-value: < 2.2e-16
##
## [1] 11
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] 12

## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.041189 -0.003543 0.000551 0.003349 0.062674
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -0.0024240 0.0001427 -16.99 <2e-16 ***
## test_set$FCH4 0.1085388 0.0078444 13.84 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.006739 on 2291 degrees of freedom
## Multiple R-squared: 0.07712, Adjusted R-squared: 0.07672
## F-statistic: 191.4 on 1 and 2291 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.088590 -0.002258 0.000217 0.001933 0.116428
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -1.279e-03 6.477e-05 -19.75 <2e-16 ***
## result$FCH4 5.009e-01 2.872e-03 174.37 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.00616 on 9169 degrees of freedom
## (18521 observations deleted due to missingness)
## Multiple R-squared: 0.7683, Adjusted R-squared: 0.7683
## F-statistic: 3.04e+04 on 1 and 9169 DF, p-value: < 2.2e-16
##
## [1] 12
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] 13

## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.047161 -0.003142 0.000342 0.003147 0.083294
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -0.0024662 0.0001453 -16.97 <2e-16 ***
## test_set$FCH4 0.0925520 0.0062235 14.87 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.006905 on 2291 degrees of freedom
## Multiple R-squared: 0.08803, Adjusted R-squared: 0.08764
## F-statistic: 221.2 on 1 and 2291 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.187353 -0.002254 0.000240 0.001963 0.105656
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -1.447e-03 7.037e-05 -20.56 <2e-16 ***
## result$FCH4 4.444e-01 3.121e-03 142.38 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.006693 on 9169 degrees of freedom
## (18521 observations deleted due to missingness)
## Multiple R-squared: 0.6886, Adjusted R-squared: 0.6885
## F-statistic: 2.027e+04 on 1 and 9169 DF, p-value: < 2.2e-16
##
## [1] 13
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] 14

## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.050861 -0.003339 0.000563 0.003200 0.084015
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -0.0023358 0.0001497 -15.61 <2e-16 ***
## test_set$FCH4 0.1272967 0.0071459 17.81 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.007116 on 2291 degrees of freedom
## Multiple R-squared: 0.1217, Adjusted R-squared: 0.1213
## F-statistic: 317.3 on 1 and 2291 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.105043 -0.002158 0.000323 0.002024 0.093241
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -1.494e-03 6.495e-05 -23.0 <2e-16 ***
## result$FCH4 4.616e-01 2.881e-03 160.3 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.006177 on 9169 degrees of freedom
## (18521 observations deleted due to missingness)
## Multiple R-squared: 0.7369, Adjusted R-squared: 0.7369
## F-statistic: 2.568e+04 on 1 and 9169 DF, p-value: < 2.2e-16
##
## [1] 14
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] 15

## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.065765 -0.003254 0.000696 0.003580 0.045389
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -0.0026641 0.0001521 -17.52 <2e-16 ***
## test_set$FCH4 0.0689746 0.0065791 10.48 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.007234 on 2291 degrees of freedom
## Multiple R-squared: 0.04578, Adjusted R-squared: 0.04536
## F-statistic: 109.9 on 1 and 2291 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.200464 -0.002251 0.000304 0.002104 0.197397
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -1.521e-03 7.842e-05 -19.4 <2e-16 ***
## result$FCH4 4.659e-01 3.478e-03 133.9 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.007458 on 9169 degrees of freedom
## (18521 observations deleted due to missingness)
## Multiple R-squared: 0.6618, Adjusted R-squared: 0.6618
## F-statistic: 1.794e+04 on 1 and 9169 DF, p-value: < 2.2e-16
##
## [1] 15
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] 16

## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.065809 -0.003324 0.000780 0.003399 0.123005
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -0.0027854 0.0001619 -17.21 <2e-16 ***
## test_set$FCH4 0.0989134 0.0070089 14.11 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.007719 on 2291 degrees of freedom
## Multiple R-squared: 0.07998, Adjusted R-squared: 0.07958
## F-statistic: 199.2 on 1 and 2291 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.146221 -0.002219 0.000435 0.002173 0.098295
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -1.684e-03 7.192e-05 -23.42 <2e-16 ***
## result$FCH4 4.527e-01 3.190e-03 141.93 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.00684 on 9169 degrees of freedom
## (18521 observations deleted due to missingness)
## Multiple R-squared: 0.6872, Adjusted R-squared: 0.6872
## F-statistic: 2.014e+04 on 1 and 9169 DF, p-value: < 2.2e-16
##
## [1] 16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] 17

## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.043002 -0.003407 0.000699 0.003223 0.038154
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -0.0028823 0.0001357 -21.24 <2e-16 ***
## test_set$FCH4 0.1220999 0.0067970 17.96 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.006467 on 2291 degrees of freedom
## Multiple R-squared: 0.1235, Adjusted R-squared: 0.1231
## F-statistic: 322.7 on 1 and 2291 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.173758 -0.002023 0.000428 0.002048 0.195287
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -1.613e-03 7.064e-05 -22.83 <2e-16 ***
## result$FCH4 5.159e-01 3.133e-03 164.69 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.006718 on 9169 degrees of freedom
## (18521 observations deleted due to missingness)
## Multiple R-squared: 0.7473, Adjusted R-squared: 0.7473
## F-statistic: 2.712e+04 on 1 and 9169 DF, p-value: < 2.2e-16
##
## [1] 17
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] 18

## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.052561 -0.003349 0.000566 0.003284 0.087575
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -0.0026149 0.0001459 -17.93 <2e-16 ***
## test_set$FCH4 0.1001950 0.0056968 17.59 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.006959 on 2291 degrees of freedom
## Multiple R-squared: 0.119, Adjusted R-squared: 0.1186
## F-statistic: 309.3 on 1 and 2291 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.198707 -0.002229 0.000508 0.002318 0.075995
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -1.757e-03 7.013e-05 -25.06 <2e-16 ***
## result$FCH4 4.119e-01 3.110e-03 132.44 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.006669 on 9169 degrees of freedom
## (18521 observations deleted due to missingness)
## Multiple R-squared: 0.6567, Adjusted R-squared: 0.6567
## F-statistic: 1.754e+04 on 1 and 9169 DF, p-value: < 2.2e-16
##
## [1] 18
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] 19

## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.053435 -0.002855 0.000555 0.002953 0.094091
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -0.0022447 0.0001495 -15.01 <2e-16 ***
## test_set$FCH4 0.1145451 0.0057200 20.02 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.007128 on 2291 degrees of freedom
## Multiple R-squared: 0.149, Adjusted R-squared: 0.1486
## F-statistic: 401 on 1 and 2291 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.175257 -0.002265 0.000385 0.002136 0.083069
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -1.582e-03 7.014e-05 -22.56 <2e-16 ***
## result$FCH4 4.157e-01 3.111e-03 133.64 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.00667 on 9169 degrees of freedom
## (18521 observations deleted due to missingness)
## Multiple R-squared: 0.6608, Adjusted R-squared: 0.6607
## F-statistic: 1.786e+04 on 1 and 9169 DF, p-value: < 2.2e-16
##
## [1] 19
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] 20

## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.042982 -0.003262 0.000772 0.003380 0.107995
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -0.0028766 0.0001485 -19.38 <2e-16 ***
## test_set$FCH4 0.1067937 0.0064476 16.56 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.00708 on 2291 degrees of freedom
## Multiple R-squared: 0.1069, Adjusted R-squared: 0.1066
## F-statistic: 274.3 on 1 and 2291 DF, p-value: < 2.2e-16




## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.190505 -0.002153 0.000478 0.002180 0.078782
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -0.0017428 0.0000693 -25.15 <2e-16 ***
## result$FCH4 0.4487030 0.0030735 145.99 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.006591 on 9169 degrees of freedom
## (18521 observations deleted due to missingness)
## Multiple R-squared: 0.6992, Adjusted R-squared: 0.6992
## F-statistic: 2.131e+04 on 1 and 9169 DF, p-value: < 2.2e-16
##
## [1] 20
## Warning: `as.tibble()` was deprecated in tibble 2.0.0.
## ℹ Please use `as_tibble()` instead.
## ℹ The signature and semantics have changed, see `?as_tibble`.
## This warning is displayed once every 8 hours.
## Call `lifecycle::last_lifecycle_warnings()` to see where this warning was
## generated.
## StartDate EndDate AvailNum EffNum ObsNum AnnNEE_gCO2
## 1 2021-09-04 2021-12-31 23:30:00 5712 535.3052 3034 464.6739
## 2 2022-01-01 2022-12-31 23:30:00 17520 1206.8784 8440 -502.5467
## 3 2023-01-01 2023-04-03 11:30:00 4440 1437.5766 2133 583.4427
## UncNEE_Rand_gCO2 UncNEE_uStar_gCO2 UncNEE_Comb_gCO2 AnnNEE_gC UncNEE_Rand_gC
## 1 42.636103 7.981962 43.37683 126.8166 11.63604
## 2 84.754166 33.961872 91.30541 -137.1526 23.13070
## 3 8.556503 8.146882 11.81463 159.2304 2.33520
## UncNEE_uStar_gC UncNEE_Comb_gC
## 1 2.178399 11.838197
## 2 9.268713 24.918638
## 3 2.223408 3.224392
## StartDate EndDate AvailNum EffNum ObsNum AnnFCH4_gCH4
## 1 2021-09-04 2021-12-31 23:30:00 5712 901.5999 965 -0.1012138
## 2 2022-01-01 2022-12-31 23:30:00 17520 4534.3640 6122 -0.2989308
## 3 2023-01-01 2023-04-03 11:30:00 4440 1201.5861 1480 -0.1774757
## UncFCH4_Rand_gCH4 UncFCH4_uStar_gCH4 UncFCH4_Comb_gCH4 AnnFCH4_gC
## 1 0.14065359 0.13619315 0.19578561 -0.07579047
## 2 0.14441551 0.07587745 0.16313561 -0.22384401
## 3 0.07083016 0.02912431 0.07658419 -0.13289658
## UncFCH4_Rand_gC UncFCH4_uStar_gC UncFCH4_Comb_gC
## 1 0.10532358 0.10198353 0.14660729
## 2 0.10814057 0.05681821 0.12215847
## 3 0.05303872 0.02180874 0.05734742